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Abstract 

The mitochondria and chloroplasts in plant cells are originated from bacteria! endosymbioses, and they 
still replicate their own genome and divide in a similar manner as their ancestors did. It is thus liltely that 
the organelle transcription is coordinated with its proliferation cycle. However, this possibility has not ex- 
tensively been explored to date, because in most plant cells there are many mitochondria and chloroplasts 
that proliferate asynchronously. It is generally believed that the gene transfer from the organellar to 
nuclear genome has enabled nuclear control of the organelle functions during the evolution of eukaryotic 
plant cells. Nevertheless, no significant relationship has been reported between the organelle transcrip- 
tome and the host cell cycle even in Chlamydomonas reinhardtii. While the organelle proliferation cycle 
is not coordinated with the cell cycle in vascular plants, in the unicellular red alga Cyanidioschyzon 
merolae that contains only one mitochondrion, one chloroplast, and one nucleus per cell, each of the 
organelles is known to proliferate at a specific phase of the cell cycle. Here, we show that the expression 
of most of the organelle genes is highly coordinated with the cell cycle phases as well as with light regimes 
in clustering analyses. In addition, a strong correlation was observed between the gene expression profiles 
in the mitochondrion and chloroplast, resulting in the identification of a network of functionally related 
genes that are co-expressed during organelle proliferation. 
Key words: cell cycle; clustering analysis; Cyanidioschyzon; organelle; transcriptome 



1 . Introduction 

The presence of organelle genomes in mitochondria 
and chloroplasts in plant cells is explained by the 
widely accepted endosymbiotic theory. According to 
this theory, mitochondria and chloroplasts originated 
from endosymbioses of a-proteobacteria-like and 



cyanobacteria-like bacteria, respectively.''^ At the 
initial stage of such symbioses, the host cell would 
have evolved some mechanism to coordinate metabol- 
ism and proliferation of the endosymbionts to control 
their growth. Subsequently, a number of genes encod- 
ing proteins for organelle functions have been trans- 
ferred to the nuclear genome, and the host cells have 
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gained the ability to dominate most organelle pro- 
cesses by controlling the proteins involved. On the 
other hand, mitochondria and chloroplasts have 
retained the ability to undergo certain proliferation 
processes such as DNA replication and division. These 
processes require coordination of events, and thus, 
these organelles must control their own transcrip- 
tomes to some degree. However, the cycle-depend- 
ence of gene expression during proliferation and the 
relevant regulatory mechanisms are not well under- 
stood for organelles, because in most eukaryotic cells 
there are numerous organelles with asynchronous 
cycles of proliferation. 

Cyanidioschyzon merolae is a thermo-acidophilic uni- 
cellular red alga isolated from an Italian volcanic hot 
spring.^ This organism has an extremely simple cell 
structure, with one nucleus, one mitochondrion, and 
one chloroplast. These three genome-containing 
organelles replicate only once at a specific phase of 
the cell cycle."^'^ Thus, this organism is an ideal experi- 
mental system to clarify the dependency of gene ex- 
pression in the mitochondrion and the chloroplast 
on the proliferation cycle. Recently, the full nucleotide 
sequences of these three genomes were deter- 
mined.^"^ Phylogenetic analyses of photosynthetic 
genes suggest that C. merolae is one of the most primi- 
tive photosynthetic eukaryotes, and that it diverged 
just after the monophyletic origin of plastids.^° The 
genome project also clarified that this alga has a min- 
imally redundant small genome. ^'^ Thus, it is a model 
eukaryote for studies on the basic architecture of 
eukaryotic cells based on interactions of orga- 
nelles."*'^ ' "^"^ Here, we report on the relationship 
between the replication cycle and the organelle tran- 
scriptomes in C. merolae. 

Most studies on transcription of the chloroplast 
genome have been performed in vascular plants. 
Their chloroplast DNA is ~150 kbp in length, and 
contains >120 genes encoding rRNAs, tRNAs, and 
proteins.* ^'^^ These genes are transcribed by the 
cyanobacteria-type plastid-encoded plastid RNA poly- 
merase (PEP). Sigma, which is the essential subunit for 
initiation of PEP transcription, is encoded by the 
nuclear genome, and often exists in multiple forms 
(e.g. six in Arabidopsis). These multiple forms are 
responsible for the PEP selectivity of promoter activa- 
tion. Also involved in this activation is a nuclear- 
encoded plastid RNA polymerase (NEP), which is 
similar to a T7-type bacteriophage RNA polymerase. 
This NEP apparently evolved from the nuclear- 
encoded mitochondrial RNA polymerase.' ^ The 
nuclear genome of vascular plants encodes multiple 
genes for NEP-type RNA polymerases, and they are 
targeted to mitochondria, chloroplasts, or both. 
Sequential actions of these two RNA polymerases 
are thought to be essential for chloroplast 



differentiation.* Genes are usually transcribed as 
operons, >60 of which have been described in 
the tobacco chloroplast genome.^* These operons 
tend to be processed and modified at a later point 
after transcription, when they are processed into 
mature transcripts competent for translation. ^^"^"^ 
The number of genes encoded in mitochondrial 
genomes differs among organisms. However, the tran- 
scriptional machinery is always composed of NEP-type 
enzymes associated with other accessory proteins.^^ 
This system appears to be rather simple compared 
with those of chloroplasts. 

Compared with higher plants, there is far less infor- 
mation available for organelle transcriptional machin- 
eries in red algae. However, it is likely that the 
organelle transcription machineries in C. merolae are 
far simpler than those in vascular plants. As revealed 
by bioinformatic analyses (http://www.cbs.dtu.dk/ 
services/TargetP/), a unique NEP-type gene encoded 
by the nuclear genome of C. merolae is predicted to 
localize in the mitochondrion. This RNA polymerase 
is likely responsible for the transcription of 62 mito- 
chondrial genes. PEP core subunits, a, (3, and (3', are 
encoded by the chloroplast genome, and in addition, 
four sigma subunits for PEP are encoded by the 
nuclear genome. Besides the RNA polymerase subu- 
nits, the chloroplast genome also encodes four 
bacteria-type transcription factors, which were likely 
inherited from the ancestral cyanobacterium. These 
transcription factors likely sustain the organelle's 
autonomous transcriptional regulation in response 
to environmental changes.^'^*^ 

In this study, we examined the relationship between 
the organelle transcriptome and the cell cycle using a 
custom-made microarray covering two organelle 
genomes of C. merolae. Previously, we constructed a 
microarray covering the chloroplast genome, and ana- 
lysed light-responsive changes of the transcriptome.^^ 
Here, we show that there is coordinated gene expres- 
sion between the chloroplast and the mitochondrion, 
and that the expressions of genes of both organelles 
are regulated both by light and by the cell cycle. In 
addition, a correlation clustering analysis revealed 
that functionally related genes show very similar 
expression patterns. 

2. Materials and methods 

2. 7. Algal material and culture conditions 

Cells of C. merolae 1 OD were cultured and their 
growth synchronized as described previously,^^ with 
minor modifications. For synchronization, cells were 
first cultured in the MA2 medium, in which the 
concentrations of all constituents with the exception 
of CaCl2, FeCia, and EDTA (disodium salt) were 
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doubled relative to the original MA medium. When 
the optical density at 7 50 nm had reached 1 0, cells 
were diluted to yield an optical density of 0.5 and 
were cultivated under a 1 8 h-dark/6 h-light 
(150|jLmol photons m"^ s"^) cycle at 42°C. 
Cultures were bubbled with 2% CO2 in air. 
Experiments with synchronized cultures were initiated 
at the end of the first dark period (1 D1 8 in 
Supplementary Fig. S2). 

2.2. RNA extraction and northern analysis 
Cyanidioschyzon merolae cells were harvested by 

centrifugation (3000 0, 4-°C, 5 min) and stored at 
-80°C until use. Total RNA was extracted as described 
previously.^" Aliquots of total RNA (5 |xg) were frac- 
tioned on a 1% agarose gel containing 0.4 M formal- 
dehyde. Probes were prepared by PCR amplification 
using specific primers designed from within each 
open reading frame. Each PCR fragment was cloned 
into the pTZblue vector (Novagen, Inc., Madison, 
Wl), and then used as a probe after digestion at a suit- 
able restriction site located in the multi-cloning site of 
the plasmid. Probes were labeled with alkali phos- 
phatase using the AlkPhos direct labeling kit (GE 
healthcare, Giles, UK) and were hybridized to RNAs 
blotted onto a positively charged nylon membrane. 
After incubation for 16 h at 55°C and suitable 
washing of the membrane, hybridization signals 
were detected using an Image Analyzer LAS-3000 
(Fuji film, Tokyo, Japan). 

2.3. DNA microarray and correlation clustering 
analysis 

We designed and constructed oligonucleotide- 
based microarrays that included all organelle 
genome-encoded open reading frames (ORFs). To 
synthesize dye-labeled probes for microarray analysis, 
we used the amino allyl RNA (aRNA) synthesis kit (ver. 
2 high yield type; Sigma). Antisense cDNA was synthe- 
sized from 5 |xg of total RNA by a reverse transcription 
reaction using a customized 55-mer primer contain- 
ing a 24-mer T7 promoter and random oligomer 
sequences (T7dT-N; SIGMA-Aldrich Japan, Komae, 
Japan). This was followed by an in vitro transcription 
reaction to incorporate amino allyl nucleotides into 
aRNA. An aliquot (1 0 jjig) of synthesized aRNA was 
coupled with Cy3 or Cy5 dyes and purified on an 
Amicon YM-30 column (Millipore, Billerica, MA). 
Probes were hybridized to the DNA microarray at 
57°C for 1 6 h, then excess probe was removed by 
washing with washing solution I (0.2 x SSC, 0.1% 
SDS) at 50°C for 1 0 min and then with washing solu- 
tion II (0.05 X SSC) at room temperature for 5 min. 

Signals on the DNA microarray were detected with 
a laser scanner GMS428 (Affymetrix, Santa Clara, 



CA) and signal intensities were calculated using 
ImaGene ver. 6.0 (BioDiscovery, El Segundo, CA). 
Signal intensity of each spot was calculated automat- 
ically with the subtraction of local background value. 
Each signal value was normalized against the averaged 
intensity of internal control spots. The microarray 
data of cycloheximide assay were further normalized 
due to relatively high local background intensity 
that causes a number of false-positive errors especial- 
ly in low-intensity signal data. Thus, we manually per- 
formed careful local background correction in each 
spot and flooring treatment by adding a constant 
value to all data of signal intensity due to Cy3 and 
Cy5 signals. The control experiment with wild-type 
cells (Supplementary Fig. SI) demonstrated that the 
range of experimental errors in the induction factor 
was between 0.5 and 2.0, as indicated by dashed 
lines in Supplementary Figs SI and S2. Clustering 
analysis of data was carried out using MeV software 
(TIGR; http://www.tm4.org/mev.html). 

2.4. Western analysis of sigma factors 

To clone genes encoding the C. merolae sigma 
factors SIG1, SIG2, SIG3, and SIG4, each open 
reading frame was amplified by PCR using C. merolae 
genomic DNA as the template. Each PCR-amplified 
fragment was inserted into pENTR™/D-TOPO 
(Invitrogen) according to the manufacturer's instruc- 
tions. To prepare the antigen for production of poly- 
clonal antibodies, each cloned gene was recombined 
into the destination vector pDEST-HIS.^° Preparation 
and purification of polyclonal antibodies and immu- 
noblot analyses were carried out as described 
previously.^°'^^ 

3. Results 

3.1 . Design of microarray for analysis of organelle 
gene expression in C. merolae 
The eukaryotic red alga C. merolae has a fully 
sequenced genome, and is an important model or- 
ganism for studying the evolution of primitive photo- 
synthetic eukaryotes.^ In this study, we monitored 
organelle gene expression in C. merolae during syn- 
chronous cultivation to understand the relationship 
between organelle gene expression and cell cycle 
events, and to determine whether gene expression 
in the chloroplast was correlated with that in the 
mitochondrion. Previously, we prepared and analysed 
a custom DNA microarray focused on chloroplast 
genes.^^ Here, we prepared a new custom DNA micro- 
array covering both the chloroplast and mitochondrial 
genomes. This enabled comparative analysis of gene 
expression in both organelles. We designed 50-mer 
oligonucleotide probes from a highly specific portion 



292 



Coordinated Organelle Transcriptome in C. merolae 



[Vol. 1 9, 



of each ORF of the chloroplast and mitochondrial 
genomes (Supplementary Table SI). These probes 
were spotted onto a glass slide as a duplicate for 
statistical reproduction. The results of a control 
experiment confirmed that experimental errors of 
our microarray system were low (within a 2-fold 
difference; Supplementary Fig. SI). 

Several groups have established synchronous 
culture systems for C. merolae.^'^^ We reconstructed 
a synchronous culture system with some modifica- 
tions that optimized cell growth in the MA2 
medium. Cultures (5 L) were grown in a thermostat 
jar fermenter and were bubbled with 2% CO2 in air 
and illuminated with white light (1 50 |xE m~^ s~^). 
Cells were first incubated for 18 h in the dark, and 
then incubated for 6 h in the light before incubation 
for a further 1 8 h in the dark. Then, cells were cul- 
tured in continuous light conditions for 24 h. The 
cell sample obtained after the first 1 8-h dark incuba- 
tion (1 Dl 8) was assigned as the time 0 point of the 
synchronization culture (Fig. 1). Cell samples were 
collected at 28 time points during the 48-h incuba- 
tion period (see upper part of Fig. 1 B and legend of 
Supplementary Fig. S2). Sample names in Fig. 1 B, 
such as 1 L0.25 (0.25 h after first light illumination), 
2D18 (1 8 h after onset of second dark period) or 
2L24 (24 h after second light illumination) corres- 
pond to 0.2 5, 24, and 48 h time points, respectively, 
in Fig. 1 A. The cellular mitotic indexes during incuba- 
tion in a light-dark cycle and under constant light are 
shown in Supplementary Fig. S2. 

3.2. Coordination of the accumulation of organelle 
transcripts with the cell cycle phase 
Previously, a genome-wide analysis reported that 
steady-state transcript levels in mitochondria do not 
vary between light and dark phases and are stable 
throughout the diurnal time course in Arabidopsis 
thaliana.^^ Another DNA microarray analysis in 
Nicotiana tabacum also showed that transcript levels 
of about half of plastid genes were not changed 
during light and dark conditions.^^ Northern analyses 
in Physcomitrella patens and in C. reinhardtii, several 
chloroplast genes showed almost constitutive expres- 
sion in diurnal conditions.^"^'^^ However, it is not 
clear that these organisms do not have synchronous 
mechanisms to regulate the transcript levels of organ- 
elle genes in a genome-wide manner or this syn- 
chrony is difficult to detect by perturbation due to 
averaged results of numerous organelles in each cell 
in individual phases of cell cycle. Since C. merolae 
has only one mitochondrion and one chloroplast, we 
anticipated that changes in organelle transcript 
levels in response to external and internal conditions 
may be more easily recognizable in this organism. 



Thus, we analysed the expression profiles of genes of 
C. merolae organelle genomes during synchronous 
cultivation (Fig. 1A). As expected, expressions of 
chloroplast and mitochondrial genes were strongly 
induced by light. We also found that a number of 
gene transcripts of both organelles showed transient 
accumulation around 4-8 h after the onset of the 
dark period. This accumulation also occurred with 
similar timing during continuous light illumination, 
which corresponded to 36-38 h after the onset of 
the first light illumination (Fig. 1 A). Thus, this expres- 
sion pattern appeared to reflect responses to internal 
events, rather than environmental light conditions. 
For both organelles, the peaks in transcript accumula- 
tion occurred at similar times. This indicated that the 
accumulation of transcripts from both organelles is 
coordinately regulated. Microscopic observation 
revealed that cells were most commonly at the 
M-Gl phases at these time points (see also 
Supplementary Fig. S2). Thus, the transient accumula- 
tions of most organelle transcripts were due to cell 
cycle-dependent regulation. 

The peak height of the mitotic index curve gradual- 
ly diminished and finally became asynchronous due 
to the diffusion of cell cycle phases among the cells 
in the culture (Supplementary Fig. S2). Similarly, ex- 
pression patterns of organelle genes also partly dif- 
fused after 42 h from the start of incubation (time 
0; Fig. 1). These results provide further evidence that 
the organelle gene expression pattern is highly corre- 
lated with cell cycle events rather than circadian 
control. In another study, qPCR analysis of the expres- 
sions of 74 selected genes (70 nuclear genes and 4 
plastid genes) during light-dark cycles'"^ showed 
that the expression of a group of light-inducible 
genes, e.g. the nuclear-encoded atpl gene for the 
alpha subunit of mitochondrial ATPase, was transient- 
ly induced after M phase during a dark period. 
Fujiwara et al.^^ also detected the accumulation of 
several transcripts of nuclear-encoded organelle 
genes, which were specific to S phase or to G2 
phase, during a dark period. These two reports indi- 
cated that phases of the nuclear cell cycle might 
affect organelle gene expression and metabolism. 
Together with our data, these results suggest that 
there is strong coordination of gene expression 
among the nucleus, the chloroplast, and the mito- 
chondrion in a cell cycle-dependent manner. 

3.3. Highly coordinated gene expression of 
chloroplast and mitochondrial genomes 
The expression patterns of chloroplast and mito- 
chondrial genes were very similar (Fig. 1 ). To estimate 
the similarity of gene expression patterns of the two 
organelles by a statistical method, we performed a 
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Figure 1. Expression profiles of organelle genes in DL and LL conditions. (A) The index of organelle division cycle (see also Supplementary 
Fig. S2) and expression profiles of 243 genes in the plastid genome and 62 genes in the mitochondrial genome are shown. The 
horizontal axis shows the duration of incubation from the onset of light illumination at starting time point of synchronization culture 
(time 0). The vertical axis shows the induction ratio of each gene compared with the basal level of expression at time 0. Expression 
profiles of chloroplast genes (upper panel) and mitochondrial genes (lower panel) were highly synchronized. (B) Hierarchical tree view 
of plastid and mitochondrial genes. Chloroplast genes (Cp), mitochondrial genes (Mt), and other functional categories are indicated. 
Each column of the heat map corresponds to each time point of sampling for DNA microarray analysis. Log2 values of the induction ratio 
of each gene compared with the basal level of expression at time 0 are shown in red (induced genes) or green (suppressed genes; 
indicator is shown at the top of the figure). Sample names in (B), such as 1 L0.25 (0.25 h after first light illumination), 2D1 8 (1 8 h after 
the onset of second dark period) or 2L24 (24 h after second light illumination) correspond to 0.25, 24, and 48 h time points from time 
0 in (A), respectively. Details of cellular mitotic indexes during incubation in a light-dark cycle and under constant light are shown in 
Supplementary Fig. S2. Functionally related genes were categorized into the same sub-clusters. The list of DNA microarray data is 
available at the following web site (http://merolae.biol.s.u-tokyo.ac.jp/download/). 



294 



Coordinated Organelle Transcriptome in C. merolae 



[Vol. 1 9, 



relative network analysis by calculating the correl- 
ation coefficient of the expression profile of each 
gene using the RelNet function of MeV software 
(http://www.tm4.org/mev.html; Supplementary Fig. 
S3). The matrix of the calculated correlation coeffi- 
cient (R^) for all genes and the list of DNA microarray 
data are available at the following web site 
(http://merolae.biol.s.u-tol<yo.ac.jp/download/). The 
results of correlation clustering showed that the co- 
expression network, > 0.7, covers 302 genes 
(98%) out of a total of 306 genes in the organelle 
genomes. This large network included almost all 
genes of organelle genomes. Expression profiles of 
only four genes, such as two rRNA genes (CMV029R, 
CMV032R) and two mitochondrial genes, 
CMW051R (trnS) and CMW002C (sdhB), were not 
categorized into this large group. Unfortunately, we 
could not identify oscillatory patterns from these 
four expression profiles. 

3.4. Highly similar expression patterns of functionally 
related genes 

In this study, we observed that genes in the same func- 
tional category showed very similarexpression patterns. 
In addition, there were only minor differences in gene 
expression patterns between the chloroplast and the 
mitochondrion. To clearly illustrate the results of the 
correlation analyses, we constructed a tree view of hier- 
archical clustering according to Pearson correlation 
coefficients (Fig. 1 B). When expression profiles of genes 
were compared, genes that were functionally related 
showed much higher correlation coefficients than 
those of unrelated genes. Thus, functionally related 
genes assembled into several sub-clusters in the hier- 
archical tree view. For instance, the tRNAclustercontains 
51 genes including 2 0 tRNAs of both organelles. The 
light-dependent gene cluster contains 2 7 plastid genes 
of photosynthetic machineries. The light- and cell 
cycle-dependent cluster contains >200 genes of both 
organelles. The housekeeping gene cluster contains 1 0 
genes such as rRNAs, RNaseP and tmRNA. The genes in 
the same sub-clusters were not from the same operon 
or gene cluster in the organelle genomes, except for a 
large gene cluster of the ribosomal proteinsand ATP syn- 
thase subunits.^ This result indicates that there is coordi- 
nated expression of functionally related genes, although 
the mechanism by which this coordination occurs is 
unknown. Mitochondrial and chloroplast genes were 
also roughly divided into distinct sub-clusters as shown 
in Fig. 1 B. 

3.5. Classification of gene expression patterns 
entrained by light or cell cycle phases 

Next, we used a template-matching approach^^'^^ 
to identify genes whose expressions were highly 



responsive to light or internal phases (Fig. 2). 
Genes were grouped by their correlation coefficients 
(i.e. the similarity of their expression to those of the 
template genes). According to the results shown in 
Fig. 1, we chose psbD and rbcL genes as template 
expression profiles for typical light-responsive and 
cell cycle-responsive genes, respectively, and con- 
firmed their expression patterns by northern ana- 
lysis (Fig. 3). Using the PTM calculation function of 
MeV software, we determined the correlation 
between the expression pattern of each gene and 
those of the two template genes. That is, we deter- 
mined whether and to what extent the expression of 
each gene resembled that of psbD or rbcL. These 
correlations are expressed as values (Fig. 2). 
The psi^D-correlated group (Fig. 2A) included many 
genes for photosynthetic machineries, and their ex- 
pression patterns were highly light dependent. 
Some of these genes consisted of operons, but 
most were located at distinct loci on the chloroplast 
genome. The expressions of the genes in this group 
were slightly induced at the end of the second dark 
period (Fig. 2A), which is likely due to circadian os- 
cillation. The rfacL-correlated group (Fig. 2B) 
included genes for translation machineries, ATP syn- 
thase, and so on. The expression of these genes 
showed transient induction just after the onset of 
light illumination as well as at the M-Gl phase 
(Fig. 2B). Some genes were expressed at relatively 
constant levels throughout the experimental period 
(Fig. 2C). This group (Fig. 2C) contained a number 
of tRNAs and several housekeeping genes, such as 
the rnpB gene. Even in this group, gene expression 
was increased or decreased at several time points 
during the cell cycle. The expression of the rnpB 
gene in Cluster III was confirmed by northern ana- 
lysis (Fig. 3). 

Next, we searched for conserved elements in the 5' 
upstream and 3' downstream regions of the genes in 
each group using the web-based tool Melina II (http:// 
melina2.hgc.jp/public/index.html). No conserved ele- 
ments were found in the 5' upstream region; however, 
we identified a conserved element in the 3' down- 
stream region of genes in each group by the 
program: CONSENSUS and MD scan in the Melina II 
[Supplementary Fig. S5 (element found in the light- 
inducible genes) and 6 (element found in the cell 
cycle-dependent genes)]. These elements were 
found in most genes in each group. The ORFs, which 
does not possess this element in their 3' downstream 
region was not the terminal gene of the putative 
operon. Positions of each element in the 120 base 
of 3' downstream region of each gene were shown 
in Supplementary Figs S5 and S6. Several genes have 
plural number of the elements in their 3' downstream 
region. 
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Figure 2. Template-matching analysis to identify light-dependent, cell cycle-dependent, or relatively constant genes in chloroplast and 
mitochondrial genomes. (A) Light-dependent genes showed a high correlation coefficient with the psbD gene, whose expression is 
highly light specific. (B) Cell cycle-dependent genes showed a high correlation coefficient with the rbcL gene, whose expression is 
highly specific to cell cycle phases. (C) Genes expressed at relatively constant levels. Upper panels show model pattern of gene 
expression used as a template for the calculation of the correlation coefficient in each group. A typical gene, which was calculated as 
highly similar pattern with that of the template pattern, was used as a second template. Lower panels show the list of genes and 
their respective values of the correlation coefficient with the expression pattern of the template genes. Middle panels show 
expression profiles of the genes, which are shown in the lower panels. The matrix of the calculated correlation coefficient (R^) for all 
genes is available at the following web site (http://merolae.biol.s.u-tol<yo.ac.jp/download/). 
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Figure 3. Northern analysis of typical genes in each cluster in Fig. 2. 
The psbD and rbslS transcripts were accumulated in a light- 
dependent manner or in cell cycle-dependent manner, 
respectively. Cells were collected at indicated time points 
during synchronous cultivation as same as Fig. 1 . The 
preparation of the probes of psbD, rbcL, and rnpB genes was 
described in Section 2. A methylene blue-stained membrane 
was shown in the bottom panel to show the amount of rRNAs. 



3.6. Expression of nuclear-encoded sigma factors of 
chloroplast RNA polymerase 
As described in Section 1 , PEP appears to be the sole 
chloroplast transcription machinery in C. merolae. 
Four nuclear genes, CMK044C, CMQ213C, 
CMR165C, and CMM072C are predicted to encode 
sigma factors for PEP, SIGl, SIG2, SIG3, and SIG4, re- 
spectively.^^ Previously, it was shown that light 
induces the accumulation of transcripts of these four 
sigma factor genes.^^ Thus, they are likely responsible 
for the light activation of chloroplast gene expression. 
To understand the link between chloroplast transcrip- 
tion and the expression of sigma factors, we analysed 
the expression of these four sigma factors at the 
protein level. First, we prepared four types of recom- 
binant proteins (SIGl, SIG2, SIG3, and SIG4) using 
an Escherichia coli system, and then obtained specific 
antibodies against each protein (see Supplementary 
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Figure 4. The expression of proteins of nuclear-encoded sigma factors for plastid RNA polymerase. Cells were collected at indicated time 
points during synchronized culture. Protein extraction and western analyses were conducted as described in Section 2. Precursors and 
mature proteins are indicated by white and black triangles, respectively. Time point 1 Dl 8 corresponds to time 0 point in Fig. 1 . The 
preparation and specificity of the antibodies were shown in Supplementary Fig. S4. The Sig2 protein showed a typical light- 
dependent accumulation. 



Fig. S4). Then, expression levels of the four sigma 
factors were determined by immunoblot analysis 
during synchronized cultivation in the same condi- 
tions as those shown in Fig. 1 (Fig. 4). In each case, 
we detected two bands, presumably the precursor 
and the processed mature protein. These bands 
were not detected in the corresponding preimmune 
sera (Supplementary Fig. S4).The size of each SIC pre- 
cursor was predicted from the C. merolae genome 
database (http://merolae.biol.s.u-tol<yo.ac.jp/), as 
follows: Sigl, 80 kDa; Sig2, 52 kDa; Sig3, 73 kDa; 
and Sig4, 67 kDa. Band sizes of CmSIGS and 
CmSIG4 proteins were higher than that predicted 
from the database. It means some modifications 
such as phosphorylation^^ of sigma factors might be 
occurring in vivo. In the case of SIG3 and SIG4, the 
band intensity of the precursor was much higher 
than that of the mature protein at several time 
points. This result suggests that sigma factor proteins 
may be regulated at the translocation and/or matur- 
ation steps. The amount of SIGl protein increased 
after the onset of light illumination and also at specif- 
ic cell cycle phases (M-Gl and Gl-S), indicating that 
SIGl may be involved in cell-cycle specific transcrip- 
tional activation. SIG2 protein was detected only in 
light conditions, suggesting that this sigma factor 
may be involved in light-dependent chloroplast gene 



expression, which is observed for genes such as psfoD 
(Fig. 3). SIG3 and SIG4 proteins were found constitu- 
tively during the dark-light (DL) and constant light 
(LL) cycles, although all four SIG gene expressions 
have been reported as light inducible.^*^ Therefore, 
expressions of SIG3 and SIG4 are probably regulated 
at translational or post-translational levels. These 
sigma factor proteins may contribute to the basal ex- 
pression of housekeeping genes in the chloroplast. 
These results indicated that the expression of sigma 
factor proteins can categorized into three groups: 
cell cycle phase-specific, light-dependent, or relatively 
constant. These three groups roughly correlate with 
the results of the clustering analysis of the DNA 
microarray. 

3.7. Light-inducible organelle gene expression 

was inhibited by cycloheximide and 

dibromothymoquinone 
Next, we examined the signaling pathways involved 
in light-inducible gene expression in organelles using 
three inhibitors: cycloheximide, 3-(3,4-dichlorophenyl)- 
1 ,1 -dimethylurea (DCMU), and dibromothymoqui- 
none (DBMIB) (Fig. 5). Cells were treated with these 
inhibitors at the end of 1 8 h dark incubation (time 
0 in Fig. 1 and Supplementary Fig. S2), and then 
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Figure 5. Effects of inhibitors on light-inducible gene expression in 
organelles. Inhibitors were added to cell cultures in dark 
conditions (1 D1 8; time 0). Cells were then incubated for 1 h 
in light conditions (1 50 yiiEm^^s"'). Total RNA was extracted 
from cells and used for DNA microarray analyses. Experiments 
were reproduced twice in biologically independent conditions 
and yielded essentially the same results. Final concentrations 
of inhibitors were as follows: cycloheximide (1 00 /xM), DCMU 
(1 0 /j,M), and DBMIB (5 /aM). Suitable concentrations of DCMU 
and DBMIB that do not affect the respiration but inhibit the 
cellular proliferation, cell division and light-induced oxygen 
evolution were determined by oxygen electrode and 
microscopic analysis. Cycloheximide and DBMIB inhibited the 
light-inducible gene expression of organelle genomes. 



subjected to light conditions for 1 fi. In an experiment 
using the cytoplasmic translation inhibitor cyclo- 
heximide, Kawazoe et a\?'^ found that light- 
inducible plastid gene expression in Chlamydomonas 
required de novo protein synthesis in the cytoplasm. 
In C. merolae, de novo synthesis of cytoplasmic 
protein was also required for light-inducible organelle 
gene expression. Taking the results shown in Figs 4 
and 5 into consideration, the inhibition of SIGl and 
SIG2 protein synthesis may be responsible for the in- 
hibitory effect of cycloheximide. We also determined 
the effects of the photosynthesis inhibitors DCMU 
and DBMIB. DCMU inhibits the electron transfer of 
Qa to Qb in the PSIl complex, inhibiting linear electron 
flow.'^" DCMU did not affect the accumulation of tran- 
scripts after the shift from dark to light. However, 
DBMIB, which is assumed to inhibit the access of 
reduced plastoquinone to the cytochrome bsf 
complex (Cyt b^f),'^^ causing the inhibition of both 
linear and cyclic electron flow, severely inhibited 
light-induced gene expression in both organelles. 

3.8. Stability of cliloroplast transcripts under ligbt 
and dark conditions 
In vascular plants, accumulation of chloroplast tran- 
scripts is usually regulated at a post-transcriptional 
level, such as stability, modification, editing, and 
translation (reviewed in).'*^ Here, we used the rifampi- 
cin assay to examine changes in stability of several 
transcripts under light and dark conditions. 
Rifampicin is a specific inhibitor of the bacterial-type 
chloroplast RNA polymerase, and has been shown to 
be effective in red algae.'^^''^'* Rifampicin was added 
to the cells 2 h after the first onset of light illumin- 
ation (1 L2) or 8 h after the second dark shift (2D8). 
As shown in Fig. 6, psbD and clpC transcripts were 
more stable during the light phase, suggesting regula- 
tion at the level of mRNA stability. However, the stabil- 
ity of rfacLS gene transcripts remained constant in light 
and dark conditions. Thus, the transient accumulation 
of rfecZ-S transcripts in the M-Gl phase (Fig. 3) during 
the dark period was likely due to transcriptional acti- 
vation. Stability of rnpB transcripts, which showed a 
relatively constant expression pattern (Fig. 3), was 
also unaffected by the light conditions. 

4. Discussion 

4. 7. Correlation between organelle gene expression 
and nuclear cell cycle phases 
In this study, we examined profiles of chloroplast 
and mitochondrial gene expression in C. merolae and 
found that gene expression in both organelles was 
highly responsive to cell cycle control and environ- 
mental light conditions. As shown in Fig. 2, genes for 
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Figure 6. Analysis of tine stability of some mRNAs in the chloroplast. Rifampicin (1 00 fig/ml) was added at indicated time point and cells 
were collected before (time 0) and after treatment (2, 5, 10, and 30 min). Signal intensity shown in lower panels was normalized 
against rRNA. Values are means of three independent experiments + standard deviation. The preparation of the probes of each gene 
was described in Section 2. A methylene blue-stained membrane was shown in the bottom panel to show the amount of rRNAs. 



transcription, translation, ATP synthesis, and RuBisCO 
are candidate genes whose expressions are obviously 
controlled by the cell cycle. However, it remains 
unclear how and why many organelle transcripts 
accumulated at the M-G1 phase, even in dark condi- 
tions. Cyanidioschyzon merolae is a photoautotrophic 
organism containing only one mitochondrion and 
one plastid, and a synchronous culture system 
entrained by light and dark cycles has been estab- 
lished. Thus, C. merolae is an ideal system for studying 
cell cycle-coordinated gene expression in organelles, 
because unlike other organelle analysis systems, it 
does not have multiple organelles with asynchronous 
replication. In vascular and other green plants, it is 
very difficult to obtain synchrony of proliferation 
cycles because of the abundant mitochondria and 
chloroplasts in the cell. In a previous study with 
tobacco BY-2 culture cells, use of organelle DNA rep- 
lication inhibitor, nalidixic acid, could synchronously 
arrest the organelle proliferation cycle, and thus 
we could successfully observe the organelle-to- 
nucleus signaling pathway and its regulatory 



mechanisms.'^''^^''^^ This kind of approach may help 
the understanding of events occurring in cells con- 
taining multiple organelles. 

4.2. Nuclear factors for light-inducible gene expression 
in the organelles 
Inhibition of cytosolic translation by cycloheximide 
abolished the light-induced transcriptional activation 
of genes from both organelles. This result suggests 
the involvement of nuclear-encoded proteins in 
transcriptional activation (Fig. 5). As described in 
Section 1 , the nuclear genome of C. merolae encodes 
only one NEP-type RNA polymerase that is predicted 
to be targeted into the mitochondrion, plastid tran- 
scription is likely dependent only on PEP Since the 
PEP activity and the specificity are dominated by 
nuclear-encoded sigma subunits, it is reasonable to 
consider that the expression of sigma factors is re- 
sponsible for the control of the plastid transcriptome. 
As shown in Fig. 4, accumulation profile of SIGl and 
SIG2 proteins showed the cell-cycle correlated 
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pattern and the light-dependent pattern, respectively. 
These patterns well correlate with the expression pro- 
files of cell-cycle correlated genes and light-depend- 
ent genes in the plastid genome, respectively, as 
shown in Fig. 2. Consistently, inhibition of light-indu- 
cible de novo synthesis of these SIC proteins by cyclo- 
heximide abolished the light-inducible expression of 
plastid genes (Fig. 5). These results strongly suggest 
the sigma factor-dependent molecular mechanism 
for the correlations provided in this study. On the 
other hand, cycloheximide did not affect the basal 
levels of transcripts in the organelles (Fig. 5). Thus, it 
is likely that the basal accumulation of transcripts is 
dependent on the constitutively expressed sigma 
factors, SIG3 and/or SIG4. Light induction of mito- 
chondrial gene expression could also be regulated by 
nuclear-encoded factors including the NEP-type 
mitochondrial RNA polymerase (CMJ2 57C) and 
other regulator(s). The effects of cycloheximide on 
mitochondrial gene expression (Fig. 5) may be 
explained by inhibition of these proteins synthesized 
de novo in the cytoplasm. 

In a green alga, C. reinhardtii which has no NEP in 
chloroplasts and it seems to have only one sigma 
factor,'^''''^^ the plastid gene expression is limited 
mainly to the post-transcriptional level. "^^'^^ Thus, 
the sigma factor in C. reinhardtii is considered as just 
a factor for basal transcription. However, our results 
suggest the red alga, C. merolae, selectively uses four 
sigma factors depending on the cellular conditions 
similarly with vascular plants. These results must con- 
tribute to the research field of sigma factor- mediated 
regulation of plastid transcription. 

4.3. Liglit-inducible organelle gene expression 
is affected by a DBMIB-related process 

What senses the light for the light-inducible organ- 
elle gene expression in C. merolae? Our results 
revealed that a process inhibited by DBMIB is essential 
for light-inducible gene expression in both organelles. 
These results, together with the results of cyclohexi- 
mide addition (Fig. 5) suggest that de novo synthesis 
of specific nuclear-encoded protein(s) is required. 
The expression of these proteins is dependent on 
light sensing around the thylakoid membrane and 
subsequent signal transduction to the nucleus. It has 
been argued that the photosynthesis-driven redox 
state is involved in light-dependent transcript accu- 
mulation in chloroplasts, i.e. the redox state is the 
mechanism for light-dependent regulation of chloro- 
plast transcription. For example, it was suggested 
that the plastoquinol/plastoquinone ratio in the 
thylakoid membrane affects the expression of the 
chloroplast psbA and psaAB genes.^' On the other 
hand, Matsuo and Obokata^^ concluded that 




Figure 7. Hypothetical scheme of coordinated gene expression in 
the plastid and mitochondrion. Nu, nucleus; Mt, 
mitochondrion; Cp, chloroplast. 

induction of the expression of Chlamydomonas PSI 
mRNAs is dependent on light illumination, but not 
on the redox state of plastoquinone. Our results indi- 
cate that the light-inducible gene expression in 
chloroplast does not depend on the redox state of 
plastoquinone pool, but depending on the integrity 
of the plastoquinol-oxidizing site of the Cytbsf 
complex. A report also suggested that a Cyt bgf- 
related factor, which is inhibited by DBMIB, affects 
the light-induced expression of several nuclear- 
encoded plastid proteins in Chlamydomonas.^^ In 
C. merolae, our results suggested that a Cyt bef-re\ated 
factor is crucial for sensing light and subsequent tran- 
scriptional activation (Fig. 7). Cyanidioschyzon merolae 
possesses neither phytochrome nor phototropin as 
photoreceptors, but has five nuclear genes that 
encode cryptochrome-related proteins. However, 
blue-light-dependent gene expression was not 
observed in C. merolae (unpublished observations), 
and thus, the photosynthetic electron transfer 
process in the thylakoid membrane may itself func- 
tion as the light sensor for transcriptional activation. 
The inhibitory effect of DBMIB is consistent with this 
scheme, whereas the signaling process from the 
chloroplast to the nucleus remains unknown. 

4.4. Correlation between changes of chloroplast 
and mitochondrial transcriptomes 
This study represents the first report of highly syn- 
chronous global gene expression profiles in the 
chloroplast and mitochondrion. This phenomenon is 
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probably due to synchronous activation of nuclear- 
encoded factors, such as SIG proteins. It is likely that 
this synchrony is essential for coordination of the 
functions of these two organelles in light- and cell 
cycle-dependent manners. The successful functioning 
of important metabolic pathways including lipid bio- 
synthesis, tetrapyrrole biosynthesis, and photorespir- 
ation rely on inter-organelle exchanges of 
metabolites. The coordinated activity of the two orga- 
nelles would be important to optimize these reac- 
tions. In another study, impaired chloroplast 
development affected steady-state transcript levels in 
mitochondria in the barley albostrians mutant.^'* 
Our synchronous system clearly showed that coordin- 
ation of gene expression between the plastid and 
mitochondrion is dynamically modulated by the 
light and cell cycle in C. merolae. 



4.5. Correlations among changes in transcript levels 
of functionally related genes 
As in bacterial transcriptional regulation, organelle- 
encoded genes tend to be transcribed as poly- 
cistronic RNAs from a limited number of 
transcriptional start sites. It is well known, especially 
in chloroplasts of green lineages, that such large tran- 
scripts are subsequently processed into smaller units, 
and that gene expression is primarily regulated at 
the post-transcriptional level. In this study, we per- 
formed cluster analysis for C. merolae organelle tran- 
scripts, and found that chloroplast genes with the 
same function frequently show similar expression pro- 
files, even when their chromosomal loci are not linked 
to each other (Fig. 2). These data indicated the pres- 
ence of a common regulator and regulatory elements 
for each cluster. As discussed above, sigma factors for 
the chloroplast RNA polymerase may function as the 
transcriptional regulators for these clusters. Recently, 
Minoda et al.^^ clarified that a plastid-encoded tran- 
scription factor, Ycf30 (CbbR), is able to bind to the 
promoter region of the rbcLS operon and up-regulates 
their expression in the isolated-chloroplasts of 
C. merolae. Together with three other transcription 
factors encoded in the plastid genome, these regula- 
tors may also be involved in the correlated gene 
expressions. In addition, we found conserved se- 
quence elements within the 3' regions of light- 
dependent (Fig. 2A) and cell cycle-dependent 
(Fig. 2B) genes. In Chlamydomonas, the endonuclease 
CSP41, which binds the stem-loop structure of the 
3' untranslated region of transcripts,^^'^^ is a major 
regulator of transcript stability. Two homologs of 
CSP41 in C. merolae, CMH047C and CML047C, are 
candidates for such post-transcriptional controls. PPR 
proteins encoded by the nuclear genome may also 
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have roles in the modulation of mRNA processing 
and the stability.^^ 

4.6. Effect of circadian rhythms on organelle gene 
expression 

Expression of photosynthetic genes, such as the 
psbD gene in the chloroplast genome, and the 
nuclear-encoded SIG2 protein showed weak induc- 
tion at the time point of artificial dawn, suggesting 
that gene expression showed a circadian rhythm 
(Figs. 3 and 4). In Chlamydomonas, nuclear regulation 
of chloroplast circadian period was observed. In 
C. merolae, it also seems that the factor due to the 
nuclear circadian rhythm affects the expression of 
SIG2 protein, and SIG2 entrain the expression 
pattern of the photosynthetic genes, such as the 
psbD with that of nuclear rhythms. In terms of cell- 
cycle control, the expression patterns of organelle 
genes appeared to be dominated by the index of cell 
division. These observations indicate that light and 
cell cycle controls are more important factors in the 
patterns of organelle gene expression than the circa- 
dian rhythm in C. merolae. 

4.7. Fundamental differences between C. merolae 
and other organisms 

In our system, synchrony of transcript levels of or- 
ganelle genes was detected in C. merolae. In vascular 
plants, such coordination of organelle gene expression 
has not been observed so far. One possible reason is 
that organelle proliferation cycle and division cycles 
are not coordinated with the cell cycle in mature 
plant cells. Second possible reason is that they do 
not have any synchronous mechanism to regulate 
the transcript levels of organelle genes in a genome- 
wide manner. Third possibility is that the synchrony 
is difficult to detect due to averaged results from 
unsynchronized organelles in each cell throughout 
the cell cycle. To discriminate these three possibilities, 
it would be helpful to check the coordination of or- 
ganelle gene expression during the specific experi- 
mental conditions where organelle proliferation 
occurs synchronously in plant cells. As Kobayashi 
etal."^ showed in tobacco BY-2 cells, organelle prolifer- 
ation synchronously begins after transferring the sta- 
tionary phased cells to new culture media. Using the 
similar experimental condition, some coordinated 
pattern of mitochondrial gene expression was 
observed by EST and microarray analyses.^" While 
only fragmental data have been obtained from vascu- 
lar plants, this kind of evidence may suggest the pres- 
ence of coordination mechanism for organelle 
transcriptome in vascular plants in specific conditions 
such as developing organ, where the mitochondrial 
DNA replication is highly activated.^' Available 
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information suggests that BY-2 cells could be a good 
model system to investigate the relationship of organ- 
elle transcription and organelle proliferation. Recently, 
Maple et al.^^ found that a nuclear-encoded gene 
whose transcript level was commonly affected in 
three arc mutants with altered plastid division pheno- 
types. This report indicates the relationship between 
plastid division cycle and cellular gene expression in 
plant cells. However, in case of mature cells of 
higher eukaryotic cells, mitochondria divide or fuse 
almost randomly and a part of them does not 
contain any DNA.^^'^"^ Thus, the relationship of mito- 
chondrial division cycle and cellular gene expression 
in plant cells is still unclear. 

Why is not organelle proliferation in other 
species more synchronous? Chloroplasts and mito- 
chondria are considered to have originated from 
bacterial symbionts and thus are possessing high 
autonomy as the nature of independent organism. 
Thus, it is not surprising that their proliferation 
and division cycles do not synchronize in a cell 
containing multiple organelles. On the other 
hand, the host cell and the endosymbionts must 
have evolved tight coordination mechanism to 
couple the organelle proliferation and division 
cycles with the cell cycle of the host cell, which 
should be most apparent when the number of or- 
ganelle is only one in a cell. Details of these coord- 
ination mechanisms are still elusive, and to be 
clarified in future studies. 
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